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LAUNCHER  DYNAMICS  STUDY 


Part  J  Phase  1 

Four-Degree-of-Freedom  Multiple  Launcher 


ABSTRACT: 

This  report  contains  an  analysis  of  a  four  degree  of  freedom  model 
multiple  launching  system.  The  launcher  pivots  about  a  point  and  the  motion 
is  resisted  by  torsional  springs.  All  of  the  missiles  are  parallel  to  one  another 
in  the  launcher  and  are  to  be  fired  one  by  one  in  a  pre-arranged  order.  The 
equations  of  motion  of  the  system  including  the  effects  of  blast  are  then  established. 
In  the  second  part,  the  effect  of  the  blast  force  of  the  missile  on  the  launcher 


face  is  considered. 
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Chapter  1 
INTRODUCTION 


1. 


This  launching  system  differs  from  the  others  because  it  contains  a  number 
of  missiles  which  can  be  fired  successively  in  a  launcher  of  box  type.  Referring 
to  Fig.  1,  we  notice  that  the  launcher  is  fixed  at  O  where  motion  is  resisted  by 
torsional  springs  and  dash  pots.  The  launching  tubes  are  parallel  and  are  open 
at  both  ends.  We  also  assume  that  there  is  no  interference  between  exhaust  cones 
of  successive  firings.  In  other  words,  at  each  time  we  consider  the  effect  of 
one  missile  only. 

The  following  information  can  be  obtained  from  this  report; 

(1)  Equations  of  motion  of  the  system  before  and  after  the  missile  leaves 
the  launcher. 

(2)  The  linear  displacement,  velocity  and  acceleration  of  the  mass  center 
of  the  missile  at  end  of  guidance. 

(3)  The  displacements  (angular)  of  the  launcher. 

(4)  An  expression  for  computing  the  magnitude  and  point  of  application  of 
the  blast  force  acting  on  a  part  of  a  ring  section. 


Chapter  II 


3. 


MATHEMATICAL  DESCRIPTION  OF  THE  PROBLEM 
1.  Coordinate  System: 


(A.B.C) 
(X.  Y.  Z) 


(x,  y,  z) 


(x_,  y_,  z_  ) 
'  m  •'m  m  ' 


(x*  ,  y'  ,  z'  ) 
'  m  m  m 


A  fixed  reference  system  on  the  earth 
Another  fixed  reference  system  on  the  earth.  Y  is  parallel 
to  the  initial  direction  of  the  launcher  tubes.  This  set 
of  axes  is  obtained  by  rotating  (A,  B,  C  )  through  a  traverse 
angle  about  the  TT-axis,  followed  by  a  rotation  of  y^ 

about  the  X-axis  (the  new  position  of  the  A-axis  after  the 
rotation  y^) 

This  set  of  axes  is  fixed  in  the  launcher  at  O  and  moves  with 
it.  The  set  is  coincident  with  (X,  Y,  Z)  initially. 

A  set  of  axes  fixed  at  the  mass  center  of  the  missile  which 
is  being  fired.  These  axes  move  with  the  missile. 

Set  of  axes  fixed  at  the  mass  center  of  an  arbitrary  unfired 
missile.  These  aixes  are  always  parallel  to  (x,  y,  z) 

Since  the  missile  always  moves  parallel  to  the  tube,  one 
coordinate  is  enough  to  describe  completely  its  motion 
relative  to  the  launcher,  r)  is  measured  from  the  mass 
center  of  the  missile  to  the  xz-plane. 


2.  Generalized  Coordinates 

This  system  has  four  degrees  of  freedom. 

The  orientation  of  the  launcher  may  be  described  by  three  successive  rotations 
^2’  ^3  as  shown  in  Fig.  1.  The  position  of  a  moving  nlissile  with  respect  to  the 


4. 


launcher  is  expressed  by  .  Thus  the  generalized  coordinates  are  0j^,  02>  0g,  ti  . 
3.  Transformation  of  Coordinates 

Following  the  same  line  of  reasoning  in  TAM  Report  139,  p.  16,  we  obtain 
the  transformation  matrices 


where 
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4.  Position  and  Velocity  Relations 

Og  :  mass  center  of  the  launcher 
From  Eqs.  (la),  (2a),  we  obtain  O  Og  in  (A,  B,  Q  as 
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Then  « 

0^  0^ 

[""2^^  ■  4" '"T”  ^  ^ 

e2  02 

+  ^b  [-"^2^1  ^  ~  ~r  "~r^~\ 

Oj^:  mass  center  of  the  firing  missile 

Following  the  same  procedure  and  taking  ±e  velocity  components  in 
(X,  Y,  Z)  we  obtain 


7. 

O  =  mass  center  of  an  arbitrary  unfired  missile. 

We  assume  here  that  i  is  the  reverse  of  the  order  of  firing,  i  lies 

between  1  and  k  .  The  (k  +  1)  missile  is  the  one  being  fired.  For  every  given 

i  ,  we  can  locate  the  missile  by  the  use  of  (x.,  ri^,  z.  ).  Note  here  ri  is  the  same 

•'  '  i  '0  i '  'o 

for  all  unfired  missiles. 


^ %  [*"2  ^^'~r  '  ^  ^  ^2  Vs 

02  q2 

+  z.  [-m2  +  ^2  ■  ~r 

i  ranges  from  1  to  k 
F  (t);  thrust  force 

We  assume  here  that  the  thrust  forces  acts  on  and  along  the  longitudinal 
axis  of  the  missile 
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The  components  of  F  are 
F„  =  -F  0a 


F  =  F 

y 


F,  =  F  0, 
z  1 


5.  Kinetic  and  Potential  Energy  of  the  System 


(10) 


Kinetic  Energy: 
T  =  Tu 


+  2T'  +  T 

m  m 


(11) 


where  T, ,  T  ,  5  T'  denote  the  kinetic  energy  of  the  launcher,  of  the 
b  m  m 

firing  missile  and  of  the  unfired  missiles  respectively. 

Tb  -1  <'ll  +  >22  ®  2  ^  '33  *3  '  '23  ^2*3  >  <"“) 

I^2»  1^3  vanish  because  of  the  symmetry  of  the  launcher 
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Potential  Energy 

V  =  V.  +  V'  +  (12) 

D  m  m  s 

where  V,.  V  ,  TV  ,  V  denote  the  potential  energy  of  the  launcher,  of 
D  m  m  s 

the  firing  missile,  of  the  unfired  missiles,  and  of  the  spring,  respectively. 
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23  2  3 


(i2d) 


Dissipation  due  to  dash  pots 

df  =  I  (£ll  +  £22  ®2^ -33  ^3  ^-12  ®1®2 

+  20^3^3+20236263) 


(13) 


6.  Equations  of  Motion 

If  0.  (i  =  1,  2,  3,  4  )  denotes  the  generalized  coordinates  (0j  ^2'  ®3  »  ’1 )» 
Lagrange's  equations  of  motion  are  : 
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The  above  equation  yields  four  differential  equations.  Three  of  these  may  be 
written  inihdicial  fottn  as  follows: 


(i,  j  =  1,  2,  3) 
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7.  Free  Motion  of  the  Launcher 

After  the  missile  leaves  the  launcher,  the  latter  will  assume  free  motion. 

Eq.  (16)  becomes  trivial  and  vanishes  identically.  The  required  equations  of  motion 
can  be  obtained  by  putting 

^  •  ^m  =  0’  ’^m  =  0*  Jil  =  Jl2  •  •  •  =  0 

^m  "  ’  ^m  "  ^T’  ^ 

The  last  term,  F  =  -T  ,  is  due  to  the  exhaust  of  the  missile  leaving  the  launcher. 

Also  note  X„,  Z_  are  quite  different  from  x  ,  z  . 

11  m  m 

The  equations  are 
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COMPUTATION  OF  THE  BLAST  FORCE 

To  calculate  the  blast  force  acting  on  the  launcher  due  to  the  exhaust  field 
of  the  rocket,  we  use  an  approximate  value  for  the  stagnation]  pressure  the 
coordinates  of  which  are  fixed  in  the  rocket.  This  force  will  be  incorrect  due  to 
the  velocity  of  the  rocket;  however,  since  we  are  also  neglecting  interaction 
between  the  exhaust  jet  and  the  launcher,  we  will  consider  it  sufficiently  accurate 
to  adopt  this  point  of  view.  The  calculation  is  based  on  the  method  outlined  in  the 
following  reference; 

"Method  of  Determining  Stagnation  Temperature-Pressure  Distribution  in 
Exhaust  Field  of  Rocket  Using  Underexpanded  Convergent-Divergent 
Nozzle"  ,  Report  Number  58-875  of  Rock  Island  Arsenal  Research  & 
Development  Division  Design  Engineering  Brapch. 

The  following  notation  will  be  used  in  this  section; 

S  the  distance  measured  from  the  nozzle  exit  of  the  missile  along 

its  axis  of  symmetry 

r  the  distance  from  a  point  to  the  axis  of  symmetry  of  the  missile 

p  (S,  r)  total  gage  pressure  at  point  (S,  r  ) 

P^  total  gage  pressure  at  the  axis  of  symmetry 

P^  pressure  at  the  nozzle  exit 


P^  rocket  motor  chamber  pressure 

y  ratio  of  specific  heats 


.17. 


M 


Mach  number  of  gas  flow  at  exit  of  nozzle. 
A. 


t 

X" 

e 


=  M 


y  +  1 


L2  +  (r-l)  M 


y  +  1 


May  be  computed  from: 


area  of  nozzle  at  throat 

area  of  nozzle  at  exit 

atmospheric  pressure 

for  S  ^  S  ,  P  =  P 
o  me 

diameter  of  nozzle  exit 
radius  of  nozzle  exit 

■"a  [<  dj  ‘'2 

-(1,  4 +  0.437  M) 
de  (S/d^ )  ^ 


(S/de  )  ^  is  given  by  Fig.  2 

the  reference  point  where  the  pressure  falls  to  l/4  of 


its  value  on  the  axis,  r^  is  given  by 


r  =  r 
0  e 


if 


S/d 


<  1 


r  =  r 
o  e 


[AtJ 
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if 


S/de 


^  1 


X  ,  z  rectangular  coordinates  in  the  plane  (surface  of  the  launcher)  normal 

s  s 

to  S  at  the  center  of  the  hole  parallel  to  the  z-axis 

outer  and  inner  radius  of  a  given  ring  section  about  the  center  of 


the  hole 


19. 


•^2 

angles  measured  with  respect  to  x  ;  they  are  the  circumferential 

s 

locations  of  the  ring  section 

T. 

1 

the  blast  force  on  the  above  ring  section 

effective  radius  of  a  ring  section 

Mx> 

M 

z 

moments  about  x  ,  z  ,  respectively,  of  the  ring  section 

2 

Erf 

error  function  =  J  e  d  X 

The  pressure  at  any  point  (S,  r)  is  given  by 


20- 


P  (S,  r)  =  a  e 


.2  2 

-b  r 


There  are  four  cases  for  the  above  expression  with  different  a,  b 


0  <S  <  S  , 
o 


r  <  ^ 


a  =  P 


b^  =  (  iiil)  2 

^e 


11 


0  <  S  <Sq  . 


k. 


a  =  P 


^2  ^  ^  2  y  2.  32  ^-2.  32 


iii 


IV 


a  =  kj  s‘'2 
e 

S  >S  ,  4^  ^  i 
0’  - 

^2 

a  =  9 

,2  .  1.  15  .  2  .2.32  _  -2.32 

b  =  (  — — )  (ka  )  (S  ) 


Using  Eq.  (18)  we  obtain 


f  fpaA.f  r 


(18) 


T.  = 

1 


P  r  d  r  d  ^ 


21. 


or 


= 


2  b 


.2„2  .2„2 

-b  R  -b  R 

(4.2-4>i)(e  '  -  e  °  ) 


(19) 


or 


R. 


<P.  jft.-  r  p  d  A 
r*fi.  r**  p  d  A 


^e- 


(  nr  r  i  ) 

.  nV  Ri>]  -  ^ V  -  Rj  ^  ) j 


1  0 
e  -  e 


d  M_  = 


■R  cos  (b 
e 


I  P  r  d  r  d  <t> 
''R; 

-  R  cos  4)  (  p  2  TT  r  d  r 


M_ 


T.  R 
i  e 


(sin  4)2  -  sin  4»i  ) 


2  TT 
2  TT 


(20) 


(21) 


Similarly 


T.  R 
1  e 


(cos  4>2  -  cos  4>^  ) 


(22) 


From  Eqs.  (19),  (21),  (22),  with  all  the  required  data,  we  can  find  the 
magnitude  and  point  of  application  of  the  force  acting  on  any  given  part  of  a  ring  section 
at  a  given  time.  Eqs.  (15),  (16),  should  yield  the  accelerations  and  velocities  of  the 
launcher  and  missile  as  the  missile  is  on  the  point  of  leaving  the  launcher.  S  may  be 
expressed  as  a  function  of  time,  i.  e.  S  -  S.  (t )  .  This  -must  be  obtained  from  the 
rocket  motion  after  firing. 
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At  a  given  time  t  ,  the  effective  area  on  the  launcher  face  covered 
by  the  exhaust  cone  (use  r  =  Tq  for  its  computation)  is  known.  We  divide  this 
area  into  a  number  of  suitably  chosen  ring  sections  with  respect  to  the  center 
of  the  hole  for  the  particular  rocket.  The  magnitude  and  point  of  application  of 
the  blast  force  on  each  ring  section  are  computed;  then  we  can  find  the  magnitude 
and  location  of  the  resultant  force.  This  will  give  T  and  X,p,  Z,^  in  Eq.  (17). 

Note  that  in  the  calculation  of  the  effective  areas,  the  effect  due  to 
holes  on  the  launcher  face  is  neglected. 
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U.  S.  Army  Artillery  Board 
Fort  Sill,  Okla. 

Department  of  Ordnance 
U.  S.  Military  Academy 
West  Point,  N.  Y. 

Commanding  Officer 
Watertown  Arsenal 
Watertown  72,  Mass. 

ATTN:  ORDBE-TX 

Commanding  General 
Army  Rocket  &  Guided 
Missile  Agency 
Redstone  Arsenal,  Ala. 

Commanding  Officer 
Chicago  Ordnance  District 
623  S.  Wabash  Ave. 

Chicago  6,  Ill. 


Commanding  General 
Aberdeen  Proving  Ground 
Aberdeen,.  Md. 

1  ATTN:  ORDBG-BRL 
1 

1  Commander 

Armed  Services  Technical 
Information  Agency 
Document  Service  Center 
Arlington  Hall  Station 

2  Arlington  12,  Vh. 

ATTN:  DSC-SE 

Commanding  Officer 
Diamond  Ordnance  Fuze 
Laboratory 

2  Connecticut  Ave.  &  Van  Ness 

1  Washington  25,  D.  C. 

ATTN:  Technical  Reference 
Section 

Commanding  Officer 
1  Rock  Island  Arsenal 
Rock  Island,  Ill. 

ATTN:  9310 
5100 

1  Commanding  Officer 
Picatinny  Arsenal 
Dover,  N.J. 

1  Director 

Marine  Corps  Development 
Center 
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a  point  and  the  motion  is  resisted  by  torsional  springs.  The  a  point  and  the  motion  is  resisted  by  torsional  springs.  The 

equations  of  motion  of  the  system  including  the  effects  of  equations  of  motion  of  the  system  includint  the  effects  of 

blast  are  then  established.  blast  are  then  established. 
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